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The G9.62+0.19-F Hot Molecular Core 

The Infrared View on Very Young Massive Stars ^ 

H. Linz, B. Stecklum, Th. Henning, P. Hofner, & B. Brandl 

Abstract 

We present the results of an extensive infrared study of the massive star- 
forming region G9. 62+0. 19. The data cover information from broad- and 
narrow-band filters in the wavelength range from 1 to 19 //m and are ob- 
tained with ESO's near- and thermal infrared camera ISAAC at the VLT and 
with the mid-infrared cameras TIMMI2 (La Silla, ESQ) and SpectroCam-10 
(Mt. Palomar). The high sensitivity and resolution provided by these facilities 
revealed intriguing new details of this star-forming region and especially about 
the embedded hot molecular core (HMC) - component F. We analyse the newly 
found infrared sub-structure of four objects in this HMC region. While one of 
these objects (F2) is probably a foreground field star, the nature of the brightest 
object in the near-infrared there (F1) remains somewhat enigmatic. Our new 
astrometry proves that this object is not coincident with the peak of the molec- 
ular line emission of the HMC, but displaced by ~ 1 .7 arcsecs which translates 
to nearly 10000 AU on a linear scale. On the basis of the available data we es- 
timate this object to be an additional embedded object with a dense dust shell. 
Very near the HMC location we find L band emission which strongly rises in flux 
towards longer wavelengths. We presume that this emission (F4) arises from the 
envelope of the HMC which is known to be associated with a molecular outflow 
roughly aligned along the line of sight. Thus, the clearing effect of this outflow 
causes strong deviations from spherical symmetry which might allow infrared 
emission from the HMC to escape through the outflow cavities. This presents 
the first direct detection of an HMC at a wavelength as short as 3.8 ^um. At 1 1 .7 
^rr\ and 18.75 jum, the HMC counterpart F4 ultimately proves to be the most 
luminous IR source within the G9.62+0.19-F region. 

In addition, within the entire G9.62+0.1 9 complex our narrow-band data and the 
K band imaging polarimetry reveal well-defined regions of enhanced Bry and 
H2 emission as well as a sector where a large contribution comes from scat- 
tered light. In combination with high-resolution radio data we make predictions 
about the extinction within this star-forming region which clarifies why some of 
the associated ultracompact Hii regions are not visible in the near-infrared. 
Our investigations show the complexity of massive star formation in full 
grandeur, but they also demonstrate that the related problems can be tackled 
by observations using the new generation of infrared cameras. 



^ Based on observations made with the ESO VLT at the Paranal Observatory under programme IDs 63.1-0329 
and 67.C-0264 and with ESO'S TIMMI2 on La Silla under programme ID 71.C-0438. 



1. Introduction 

During the last 10 to 15 years, a canonical picture of the formation of low-mass stars has been estab- 
lished (Shu et al. I1987L Andre et al. |T993|), although details of the formation process are still under 
discussion (e.g., Hartmann .20031) . On the other hand, a corresponding consistent model for the birth 
of massive stars (M > 8-10 M©) has still to be developed (see, e.g, reviews by Kurtz et al. 20001 
Stahler et al. I2000L Yorke I2003D . On average, sites of high-mass star formation are found at larger 
distances from the Sun than regions of low-mass star formation. In addition, they are located in more 
crowded and more heavily obscured regions than young low-mass stars. Both facts strongly hamper 
observational investigations. Furthermore, the evolution from a contracting cloud core to a high-mass 
main sequence star occurs on much shorter time scales, whereby the protostar shows more intensive 
interaction with the natal environment which further increases the complexity of the problem. Thus, 
the requirements to study massive star formation are severe for both, the observer and the theoretician 
(Garay & Lizano fT^^ Henning et al. 2002). 

Over the last decade, it has become possible to explore earlier stages in the evolution of massive 
young stellar objects (YSOs). While for many years ultracompact HII regions (UCHlIs) played the 
central role in these investigations (e.g.. Wood & Churchwell 1989, Kurtz et al. 119941 ). recent stud- 
ies have revealed that warm (T > 100 K), dense (uh, = 10^..10^cm^'') and compact (size < 0.1 pc) 
condensations within the molecular clouds, the so-called hot molecular cores (HMCs), can provide 
insights into even earlier phases of the formation of high-mass stars (Cesaroni et al. 119941 Kurtz et 

ai.Eunni). 

It was found that HMCs are close to, but often not coincident with, the adjacent UCHlIs (Cesaroni 
et al. 119941) . Several authors have proposed an evolutionary scenario leading from HMCs to UCHlIs 
and, finally, to more evolved HII regions (Cesaroni et al. II 9941 Kurtz et al. 120001) . 

The well-studied galactic region G9.62-l-0.19 (Garay et al. ll993[ Hofner et al. 119961 ) at a distance 
of ~ 5.7 kpc (Hofner et al. 1994) contains a number of massive YSOs in different evolutionary stages 
located close to each other. In this respect, there are indications that an age gradient goes from western 
(older) regions to eastern (younger) ones (Hofner et al. 1 19941 [19961 Testi et al.fT^8). One of the most 
interesting objects of the region (source F) is confirmed to be an HMC (Cesaroni et al. 1994, Hofner 
et al. 1 19941 1 1996") by means of interferometric observations of the dense gas tracers NH3 and CH3CN. 
Situated between two UCHlIs (radio components D and E), weak cm-continuum radiation was found 
recently for component F using high-sensitivity VLA measurements (Testi et al. 2000) . 

Another feature, which makes this HMC a unique target, is the fact that very close to the loca- 
tion of G9. 62-1-0. 19-F emission in the 2.2 /urn continuum (K band) was detected (Testi et al. 119981) . 
Corresponding to the standard models of HMCs, describing them as spherically symmetric and ex- 
tremely dense dust and gas spheres (Kaufman et al. 1 998l Osorio et al. 11999 1 the optical depth 
should be so high that the HMC should remain undetectable at near- and mid-infrared wavelengths. 
This raises the question, what the actual origin of the K band emission is. In this respect it is impor- 
tant to mention that Hofner et al. (120011) have detected a molecular outflow in G9. 62-1-0. 19 with the 
flow orientation very close to the line of sight, for which they favour component F as the most likely 
driving source. The working hypothesis of Hofner et al. (120011) was that the outflow forms a cavity of 
strongly decreased density roughly along the line of sight through which the infrared (IR) emission 
from the HMC can escape. 

While the G9. 62-1-0. 19 region was covered over the years by many radio observations with more and 
more increased resolution and sensitivity, detailed information about this star-forming region in the 
infrared is rather scarce. By utilising high-resolution near- and mid-infrared data, we have explored 



the IR properties of this region. In Sect. 2 we explain details of the observations we have conducted 
during the last four years. In Sect. 3 we present the direct results of these measurements and discuss the 
global features of the region and the specifics of the HMC zone on which the emphasis of our analysis 
will lie. Since an accurate astrometry will be a crucial issue for the following interpretation of the data, 
we included a detailed description of the astrometric techniques we have applied. In combination with 
interferometric radio maps and a radiative transfer model, we discuss the implications of our infrared 
data for an interpretation of the objects in the HMC region in Sect. 4. Furthermore, we compare the 
HMC region to the UCHlIs nearby and to the HMC in W3(H20). Finally, in Sect. 5 we summarise 
our conclusions. 

2. Infrared observations 

Table 1. Summary of observations 



Date 


Filter and A^ 


Telescope / Camera 




FOV^ 


PSF FWHM 


Ref. Star 




1999 Sep 


J 


1.25 jum 


ANTU / ISAAC 




151" 


0'.'55 


S234-E, S071- 


-D 


1999 Sep 


H 


1.65 yum 






151" 


a'65 


S234-E, S071- 


-D 


1999 Sep 


Ks 


2.16 jum 






151" 


a'55 


S234-E, S071- 


-D 


2001 May 


nb_K 


2.09 jum 






M51" 


^'65 






2001 May 


H2(1-0)S1 


2.122 jum 






151" 


0'.'48 






2001 May 


Bry 


2.166/mi 






151" 


a'44 






2001 Jun 


V 


3.78 jum 






79" 


^'55 X ^'42 


HR 6070 




2001 Jun 


Bra 


4.07 jum 






79" 


a'48 X ^'41 


HR 6070 




2001 Jun 


nb_M 


4.66 jum 






79" 


0'.'47 X ^'40 


HR 6070 




2003 Jul 


Nl 


8.70 jum 


ES0 3.6-m/TIMMI2 


64' 


x48" 


0'.'86 


HD 169916 




1999 Jun 


N 


11.70jum 


5-mHale/SC-lO 




16" 


^'76 


5 0ph 




2003 Jul 


Q2 


18.75 jum 


ES0 3.6-m/TIMMI2 


64' 


x48" 


1"60 


HD 187642 





^ Field of view for single frames, respectively 

'' Imaging polarimetry with Wollaston prism and slit mask in the optical path 



2.1. J, H, and Ks band observations 

The near-infrared data were obtained in September 1 999 with the infrared camera ISAAC (Moorwood 
1 1 9971 1 at the VLT ANTL/ telescope within the programme ID 63.1-0329. We used the opportunity to 
utilise these archived data for our purposes. 

The following NIR broad-band filters were used: J {A,- = 1.25//m, FWHM = 0.29 //m), H {Ac = 
1.65 yum, FWHM = 0.30 yum), K^ (A^ = 2.16//m, FWHM = 0.27 //m). The seeing conditions were 
between 0.65 arcsec (H-band) and 0.55 arcsec (J- and Kg-band), respectively. The pixel scale is 
O'.'147/px which results in a field of approximately 2'.5 x 2'.5. Within the overlap area of the 5 offset 
images per wavelength, the total integration time on source for each filter was 266 s . After the stan- 
dard procedure (dark-field and flat-field reduction, bad pixel removal and sky subtraction), utilising 
our own IDL^-based reduction pipeline for the ISAAC data (see also Stecklum et al. 120031 and Apai 

^ Interactive Data Language 



et al. '2004'), we corrected for the field distortion by applying the correction terms recently provided 
by ESO. Then the 5 dithered frames were combined to a larger mosaic image achieving subpixel 
accuracy. As standards for photometric calibration, the stars S234-E and S071-D^ (from the list of 
faint NIR standard stars by Persson et al. 1998| were observed before and after the observation of 
G9. 62+0. 19, respectively. 



2.2. Br y and H2 observation and K narrow-band polarimetry 

Images of G9.62+0.19 in the Br y (A^ = 2.16//m, FWHM = 0.02 /urn) and in the H2(1-)S1 (Ac = 
2.12/zm, FWHM = 0.02 /um) narrow-band filters were obtained in May, 2001, again with ISAAC 
at the VLT ANTU telescope. Pixel scale and field of view are the same as for the JHKs broadband 
imaging. The total integration times on source for both filters were 500 s within the overlap area of 
the 5 off"set images per wavelength. 

The imaging polarimetry system of ISAAC was used to obtain complementary NIR data for the 
G9. 62+0. 19 region. The camera is equipped with a WoUaston prism which divides the infalling light 
into two perpendicular polarised beams having a separation of ~ 22". The prism is used with a special 
slot mask (3 opaque stripes, each being 20" wide) to keep the two fully polarised copies of the imaged 
area separated on the detector. The first step is to simultaneously obtain the 0° and 90° data. Then, 
the reference plane is rotated by 45° by offsetting the image derotator of the alt-az mounting control 
system of the telescope accordingly. In this way, the 45° and 135° data can be obtained. To ensure a 
seamless coverage of the entire star-forming region, we used a 3-point dither pattern perpendicular 
to the stripe orientation. A narrow-band K filter was used (Ac = 2.09 yum, FWHM = 0.02 yum), which 
does not contain strong spectral lines. 



2.3. L and M band observations 

ISAAC is also capable of imaging in the U (A^ = 3.78 jum, FWHM = 0.58 jum) and nbJM 
(Ac = 4.66 /im, FWHM = 0.10 //m) thermal infrared narrow bands. Within our programme 67.C- 
0264, service time observations were performed for 09. 62+0. 19 in June, 2001, thus, very shortly 
after the repair of the ISAAC chopping secondary mirror in April and May, 2001 . The fact that during 
the first weeks after recommission of the chopping mode the data were taken without field stabilisa- 
tion resulted in a slightly degraded image quality, recognizable by means of slightly elongated star 
shapes. In addition to the above-mentioned filters the 4.07 //m narrow-band filter (FWHM = 0.07 
jum) was used which includes the Brackett a line. The pixel scale in connection with the ALADDIN 
detector array is O'.'071/px, the field of view is then 73" x 73". Since the L' and nbJM bands lie within 
the thermal infrared, on-source/off-source chopping is required to remove the thermal background 
in the images. Therefore, we have the positive and the negative beam in the images, separated by 15 
arcsec. In addition, a 5 point dither pattern was applied. The standard star HR 6070 (van der Bliek et 
al. 119961) was used to calibrate the data. 



^ We mention here that, unlike the other finding charts shown in Persson et al. (I1998t . the finding chart for 
S071-D is (accidentally) not centred on the actual calibration star which at first caused some confusion for our 
calibration. 



2 A. N and Q band observations 

We utilised two different cameras to image the central region of G9.62+0.19 in the mid-infrared 
(MIR). The 1 1.7 /im image was obtained in June 1999 using SpectroCam- 10 (Hay ward et al. 1993) at 
the 200-inch Hale Telescope of the Palomar Observatory"^. The camera uses a 128" pixel Si: As BIB IB 
detector manufactored by Rockwell. For imaging the chip is binned into a 64 x 64 pixel array with 
O'.'256/pixel. The effective wavelength of the filter was 11.7 //m with a FWHM of 1 //m. Chopping 
and nodding was performed in the standard way (chopping parallel to nodding) with a throw of 20 
arcsec. After applying a standard chop-and-nod reduction for the raw frames, a 5-frame mosaic was 
combined into the final image. The average on-source integration time at each pixel is 200 s. The star 
S Oph (HD 146051) served as standard for the flux calibration. 

Second, we used the TIMMI2 camera (Reimann et al. 120001) at the ESO 3.6-m telescope for additional 
MIR observations in July 2003. Equipped with a 320 x 240 px^ Si: As BIB Raytheon array, the FOV 
in the high-resolution mode is 64" x 48". Hereby the Nl filter {A^ = 8.7 //m, FWHM = 1 //m) and 
the Q2 filter (Ac = 18.75//m, FWHM = 0.9 yum) have been used. The chopping and nodding throws 
were 15 arcsec, but both movements are performed perpendicular when using the TIMMI2 standard 
setup. After combination of the chopped and nodded beams the total on-source times are 749 s for 
Nl and 795 s for Q2, respectively. While the sky was stable and almost photometric during the Nl 
measurements, the Q2 band observations on the following day suffered from mediocre atmospheric 
conditions. The stars HD 169916 (Nl) and HD 187642 (Q2) from the list of TIMMI2 standard stars 
were used for calibration. We mention that the Q2 standard star was observed immediately after the 
end of the Q2 exposure for the science target und hence was affected by the same sky conditions. 

With the TIMMI2 data we also tried out a different (and perhaps more sophisticated) approach 
to create the final image, namely the so-called projected Landweber matrix restoration (Bertero et 
al. 1^000 and references therein). We refer to appendix Id for some comments about the restoration 
method itself and the modifications we applied. 

Because of our modified chop/nod observational approach we have to deal with artifacts (see ap- 
pendix) in two orthogonal directions. We have undertaken the first steps to reduce these collateral 
effects (e.g., ghost images) by precisely aligning the array orientation along the chopping/nodding 
throws and by carefully registering the real values of the throws. A simple way to minimize the mul- 
tiple artifacts would be, according to Bertero et al. (120001) . to combine several restored images each 
having different throws. Due to time constraints, however, we could only perform one long integration 
for each filter, so we were not able to completely avoid these perturbing effects. In this respect one 
has to remember that a part of the original signal is then still contained in the artefacts. The effects 
on photometric accuracy are under investigation. For this paper we will use the restored TIMMI2 
images only to reveal the MIR morphology of the region. The photometry was always performed on 
unrestored and undeconvolved images. 

3. Results 

3.1. NIR data 

In the Figures IC.ll and IC.2I we show 3-colour-composites derived from the VLT J-H-Kg data, 
whereby Figure IC.2I shows the inner part of Figure IC.ll in more detail and serves as a reference 
for the radio observation nomenclature we have adopted to indicate the various components. 

^ Observations at the Palomar Observatory were made as part of a continuing collaborative agreement between 
the California Institute of Technology, Cornell University, and the Jet Propulsion Laboratory. 



Figure EU displays the global morphology of the star-forming complex G9. 62+0. 19. The extended 
diffuse region in the centre of the image dominates the near-infrared emission and can be related to 
the more evolved HII component B in Fig. IC.2I This region harbours a cluster of modestly embed- 
ded stars. Many of them show infrared excess. Recently, Persi et al. (2003) analysed in detail this 
population of stars in the G9. 62+0. 19 B and C regions. 

Eastwards of this emission, one can notice an extinction gradient. Hence, the neighbouring area, 
beginning a few arcsec away from the eastern border of component B, appears darker and more 
obscured. The superimposed MSX contour lines which we derived from the 8.28 yum image of the 
related MSX source (Egan et al. 119961) experience the same decrease as the NIR emission towards 
the eastern darker region. This is a hint for a similar behaviour of the optical depth in the near- and 
mid-infrared. In fact, we categorise this eastern region as a so-called Infrared Dark Cloud (IRDC - 
see Egan et al. .l998j . 

The transition stripe between the two above-mentioned regions is aligned roughly along the north- 
south direction and is populated with several more compact radio sources named C, D, E and F 
(Hofner et al. 119961 see also Fig. lC.2l for reference). 

Mainly in the southwest of the radio component C we detect strong K band emission. The peak 
positions of the 3.6 cm and the NIR emission are clearly separated by almost 3 arcsec. This shift is 
obvious in Fig. lC.21 it is real and not an effect of an unprecise astrometry. The K narrow-band data 
and polarimetry data (Secs. l34l and l3.5l) will enlighten these circumstances. 

The components D and E are ultracompact HII regions (Hofner et al. 1996" Testi et al. 12000") . 
whereby E might be associated either with an expanding ionised shell or with a wide-angle molecular 
outflow (Minier et al. i2001i) . Especially source D has a strong peak at 3.6 cm which can be recognized 
via the superimposed VLA contours in Fig. lC.21 Both components show no counterparts in our NIR 
data. We will discuss the consequences of these findings in Sect. 14.31 

The HMC component F was originally defined as the peak in NH3(5,5) and NH3(4,4) transition 
measurements, conducted with the VLA (Hofner et al. 1 1 9941 Cesaroni et al. 1994). The black plus 
sign in Fig. IC. 21 pinpoints the peak position of this molecular line emission. (In Fig. IC.3l we over- 
lay the contours of the HMC NH3(5,5) emission onto the NIR data for comparison.) Further radio- 
interferometric observations could reveal emission peaks in the 2.7 mm dust continuum (Hofner et al. 
119961) as well as in the cm free-free continuum (Testi et al. 120001) at the position of the HMC com- 
ponent F Testi et al. (119981) have reported a K-band detection at the coordinates of the HMC. Our 
infrared imaging reveals an intriguing structure in this HMC region which one can decompose into at 
least three distinct objects. We refer to Fig. lC.41 for the numbering we have chosen for the objects and 
which will be used in this paper. On the basis of our high-resolution VLT data (pixel scale = (X.'lS/px, 
seeing Qf.'SS) one can presume that this previous K-band observation of Testi (pixel scale = O'.'92/px) 
was suffering from the lower resolution and confusion among the objects in that line of sight which 
can now clearly be separated. We refer to our appendix lAl where we discuss in detail our approaches 
to derive an accurate astrometry for our data. 

Near the HMC peak coordinates, we find a "red" compact source according to the 3-colour-composite 
(object Fl in Fig. lC.41) . However, Fl does not coincide with the hot core peak F, but is displaced by 1 .7 
arcsecs. Although this finding suggests that the compact NIR emission Fl is not a direct trace of the 
innermost HMC region but a distinct object we should note that according to Cesaroni et al. ( 1994 ) the 
warm ammonia emission is extended over several arcsec and a source size of 3 - 5 arcsec was derived. 
Since 1 arcsec translates to 0.028 pc (assuming a distance of 5.7 kpc) Fl might still be affected by 
warm and dense molecular gas of the actual HMC. Thus, we have derived the Kg band magnitude of 
object Fl by using S234-E (K^ = 12.070 mag) and S071-D (K^ = 11.839 mag) (Persson et al. lT9981) 



as photometric standards. By performing a PSF photometry on the mosaiced science frame we derived 
a count rate for our object. To be in accordance with the 10" aperture photometry applied by Persson 
et al. (IT998 ) we also used this method (in form of the MIDAS^ magnitude/circle routine) to derive 
the count rates for the standard stars. Corresponding to this procedure the Fl magnitude is Kg = 13.6 
mag with 0.1 mag as a conservative error limit. Testi et al. ( 1998) had derived K= 12.9 ± 0.08 mag, 
probably for all the 3 objects visible at 2.16 fim in the F region in Fig. lC.4l 

At the position of the HMC 3.6 cm radio peak we find another much fainter compact source (object 
F2 in Fig. lC.4l) . In order to evaluate this fact one should remember that this area near the IRDC is in 
general very extincted. It is possible to derive a preliminary extinction estimation for this region from 
the data contained in Schlegel et al. (119981) which results in A^ values between 31.2 mag (J band) 
and 12.7 mag (K band). Therefore, other objects in the F region appear to be strongly reddended, 
with (H-K) values clearly exceeding 1.5 mag. However, this object has a "yellow" colour index with 
(H-K) being only 0.72 mag. By inserting the data for F2 into a JHKs two-colour diagram we see that 
the F2 colours can be explained almost completely by interstellar reddening. Therefore, we conclude 
that F2 is a foreground star, not physically related to the star-forming region. A definitive answer 
cannot be given until high-resolution NIR spectroscopy is performed with an instrument capable of 
resolving the objects in the HMC region (e.g., NACO at the VLT). 

A third very faint source (object F3) is visible less than 1" to the north of the foreground star. This 
object appears very red because it is only detected in the Ks band and not in the J or H band. An 
interesting fact to mention is that, according to our astrometry, an H2O maser (the first one in the H2O 
maser list in Tab. El) is only 0.25 arcsec apart from this object's position. 
The photometry of the objects in the HMC region is summarised in Table El 



3.2. L and nbM band data 

This wavelength range is interesting because it fills the gap between the Kg and the N band. Naturally, 
the signals of most of the NIR-blue field stars strongly fade and finally disappear. Hence, they do not 
cause confusion anymore. Still, enough objects remain in the FOV of ISAAC to perform an accurate 
astrometry. Hereby we want to direct the attention to an interesting detail (see Fig. lC.4l) . 
In the Ks band image one can clearly see the presence of three objects in the HMC region F. The "red" 
object Fl, already mentioned in subsection 13. 11 is dominating. The object F2 which we estimated to 
be an unrelated foreground star is almost 1 .7 arcsec apart in the north-east direction. Above of this 
we detect the third very red and weak source F3. Please note that the location of the foreground star 
is underlaid with weak and diffuse emission. 

When we now switch to the L' band we notice that the point-like star F2 in the middle has almost 
totally disappeared. Instead, we now clearly see emission arising from a slightly different location, 
that reveals the presence of another object F4. This is evident, when we compare the relative positions 
of the emission centres. We mentioned the distance of 1.7 arcsec between the point-like star F2 and 
the "red" object Fl in the Kg band. The centre of F4 in the L' band image is clearly shifted compared 
to the F2 position; it is now only 1.3 arcsec away from object Fl. This shift is persistent also in the 
4.07 jum and the nb31 band image; it is a small effect, but noticeable. In the following, the nature of 
object F4 will become more clear. 



Munich Image Data Analysis System 



Table 2. Compilation of positions for features in the HMC region 



Feature 



RA(J2000) DEC(J2000) Telescope 



NH3(4,4) 


18 06 14.85 


-20 3138.91 


VLA - D(') 


NH3(4,4) 


18 06 14.81 


-20 3139.41 


VLA - C^i) 


NH3(5,5) 


18 06 14.87 


-20 3139.41 


VLA-C(2) 


CH3CN(6-5) 


18 06 14.89 


-20 3139.20 


0VR0(3) 


2.7 mm cont. 


18 06 14.89 


-20 3138.90 


0VR0(3) 


3.6 cm cont. 


18 06 14.88 


-20 3139.37 


VLA - BnA^'^) 


H2O masers 


18 06 14.88 


-20 3138.39 


VLA-B(5^ 




18 06 14.87 


-20 3141.21 






18 06 14.81 


-20 3137.43 






18 06 14.93 


-20 3141.55 






18 06 14.86 


-20 3140.59 






18 06 14.78 


-20 3134.71 






18 06 14.87 


-20 3137.57 






18 06 14.67 


-20 3131.44 






18 06 14.98 


-20 3142.53 




object Fl 


18 06 14.78 


-20 3140.48 


ISAAC^^^ 


object F2 (Ks) 


18 06 14.87 


-20 3139.50 


ISAAC^^^ 


object F3 


18 06 14.87 


-20 3138.39 


ISAAC^^^ 


object F4 (nb_M) 


18 06 14.85 


-20 3139.70 


ISAAC^^^ 



Ref.: (1) Cesaroni et al. (fT994l) , (2) Hofner et al. (fT994l , 
(3) Hofner et al. (fT996b , (4) Testi et al. (EnHHl , 
(5) Hofner & Churchwell (fT996l , (6) this article 



3.3. N and Q band data 

Fig. IC.5I shows the central part of Fig. IC. II but now at A = 11. 7 /zm obtained with SpectroCam-10. 
To amplify the morphological structure, we applied 25 iterations of a modified Maximum Entropy 
Method for noise suppression, based on the wavelet transform and developed by Pantin & Starck 
(119961) . The extended emission, coming from component B, dominates also in the mid-infrared. A 
second, but compact source is found near the extended B component. The accurate position and hence 
the correct interpretation of this compact source has been a matter of debate during recent years (cf . 
De Buizer et al. 120001110011 Stecklum et al. lMTTl Persi et al. l2003l) . In AppendixEwe describe the 
circumstances and explain the result of our new astrometry. We find that the compact MIR emission 
is coming from the HMC region, contrary to the results of De Buizer et al. ( 2003 ). We identify it with 
object F4 (see Fig. lC.4l) . i.e., with the HMC conterpart.^ For a determination of the flux coming from 
the compact source F4, we observed the star S Oph just before we imaged the G9. 62-1-0. 19 region. 
This isolated star is a strong IRAS source with a 12 yum flux of 149.7 Jy. By applying an aperture 
photometry to the compact source (using a 3 arcsec diaphragm) and the photometric reference star, 
we derived a flux of (0.60 + 0. 10) Jy for F4. 



^ Persi et al. C2003^ also identified the MIR source with the HMC on the basis of low signal-to-noise data 
and related it to the K band emission found by Testi et al. ( 1998 1. However, they apparently underestimated the 
significance of the multiplicity of objects in the G9.62-i-0.19-F region and thus could not distinguish between 
Fl (dominating in the K band) and F4 (dominating in the N band). 
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The 8.7 //m Nl data proved their value especially regarding the MIR astrometry, as shown in Appendix 
IbI We show these data in Fig. lC.6l To produce this image the matrix restoration method described in 
appendix O was applied to the Nl band TIMMI2 data. Thereafter, the restored data were overlaid as 
contours on the NIR colour composite. As can be easily seen, they comprise a much larger field of 
view than the SpectroCam-10 data in Fig. IC.5[ We also derived the flux of the compact MIR emis- 
sion at this wavelength of 8.7 //m. As is evident in Fig. IC. 61 and Fig. IC. 41 the emission is elongated; 
most of the flux stems from object F4, but object Fl might still contribute a small fraction to the total 
signal. We did not try to disentangle F4 and Fl and only measured the total flux with an adequately 
sized aperture for the photometry. One should keep in mind that the Nl filter covers a prominent 
R\H feature, the related emission also peaks around 8.6 - 8.7 /zm (e.g., Peeters et al. l2002bl) . In fact, 
such R\H emission seems to be a common feature of the majority of compact H II regions and their 
surroundings (Peeters et al. 2002a) . ^ Thus, the measured Nl flux is probably affected by additional 
contributions from those small PAHs that are suscptible to transient heating and often attain a higher 
(non-equilibrium) temperature than the larger normal dust grains. Indeed, the Nl flux of (1 . 10 + 0. 10) 
Jy is higher than the 11.7 //m flux of F4. 

Although the atmospheric conditions during the Q2 band observations at 18.75 yum were not photo- 
metric we tentatively report the measured flux of (22.8 + 7.0) Jy for the compact emission blob we 
assume to be identical to F4 - no other reference objects appeared in the data except of the extended 
emission of component B. The quoted relatively high uncertainty of 30 % is mainly caused by the 
non-perfect sky. 

All the MIR measurements are also included in the SED plot of object F4 in Fig. IC. 91 and draw an 
enlightening picture of the energetics of the HMC. In Sect. 14.1.11 we will further elaborate on the 
consequences of our findings for object F4. 

3.4. H2 and Bry narrow band imaging data 

In Fig lC.7l we combined the Br y data (blue channel) and the H2 data (red channel). The average of the 
signal of both filters was put into the green channel. The Br y emission of the compact Hii region B 
clearly dominates the overall appearance of the star-forming complex. The ultra- and hypercompact 
Hii regions D and E are not visible at all due to the large K band extinction towards these components 
(see section |431 for a discussion). More interesting, the signal arising from the region of the radio C 
component appears almost white in this colour composite; thus, it represents just the contribution of 
the K band continuum without strong H2 or Br y line emission. The same we find for the NIR objects 
around the HMC region F 

The H2 data reveal a compact emission region located between the components B and 
C. As a whole, it appears as a flat ellipsoid with a large axis of about 5", but it is dis- 
sected into several emission maxima and minima. This feature was also mentioned by Persi 
et al. (20031) and we report here the positions of the two main components according to 
our astrometry.*^ Component 1: ajooo = 18''06™14H2,52ooo = -20°31'36'.'8, Component 2: 
0^2000 = 1 8*^06" 14';34, (52000 = -20°31'36'.'8. The projected distance to the hot core region F is 
roughly 7 arcsec (i.e. 0.2 pc if assuming d = 5.7 kpc). In principle, such a distance would allow 
an association with the molecular outflow arising from component F (Hofner et al. .2001:) . One 



^ We mention our experience with another massive star-forming region (Apai et al. 12004 !!) where PAH emission 
(also covered by the MSX A band) clearly affects the SED of the central object. 

^ For the probable cause of the astrometric disagreement between our positions and the one of Persi et al. 
(120031 we refer to the end of Appendix 1X1 
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Table 3. Measured IR fluxes in Jy (corresponding magnitudes in parenthesis) for the objects in the 
HMC region 



Band name 


object 


object 


object 




object 




Fl 


F2 


F3 




F4 


J 


- 


2.21 xlO-"* (17.15) 


- 




- 


H 


6.19x10-5(18.00) 


5.25 x 10-4 (15.72) 


- 




- 


Ks 


2.25 X 10"^^ (13.60) 


6.57 X 10-4 (15.00) 


1.67 X 10-4 (16.49) 




- 


U 


2.43 X 10-2 (10.04) 


- 


1.00x10-3(13.50) 


2.95 


X 10-3 (12.33) 


nb_M 


3.10x10-2(9.21) 


- 


- 


1.49 


X 10-2 (10.01) 


Nl 


1.10 ( 3.90)^ 


- 


- 


1.10(3.90)^ 


N 


- 


- 


- 


6.00 


xlO-i (4.55) 


Q2 


- 


- 


- 


2.28 


X lOi (-0.85)'^ 



^ In the undeconvolved Nl band, the objects Fl and F4 merge to one entity, whose combined flux is reported. 

^ Note that the Q2 flux is affected by a relatively large uncertainty of around 30% due to non-photometric sky 
conditions (see text). 



problem is that the outflow is well aligned along the line-of-sight. Thus, if the interpretation of 
shock excitation due to interaction with the outflow holds, then the outflow driving source and the 
H2 emission region cannot be located in the same plane of the sky, implying that the latter might 
be detached from the actual star-forming complex G9. 62+0. 19. Another possibility is the existence 
of another outflow in this region that was not detected in our previous studies (Hofner et al. 119961 
1^00 1 ). An alternative consideration would be that a considerable fraction (if not all) of the detected 
H2 emission is in fact fluorescence, excited by UV photons of the nearby Hii regions (e.g.. Black 
& van Dishoeck 119871 and Draine & Bertoldi [19961 for a theoretical treatment, e.g., Fernandes et al. 
I1995l and McCartney et al. 1999 for observational evidence in the case of other astronomical objects). 
The components B or C could act as UV photon donor. Only further spectral analysis, comprising 
several H2 roto-vibrational transitions, can distinguish between the coUisional shock model and the 
fluorescence model. However, we mention that the H2 feature is not a pure emission line object, we 
also see diff"use continuum emission from this position in the L' and nbJM band, as well as in the 
2.09 [xm filter used for the polarimetric imaging (see the directly following paragraph). 



3.5. K narrow band imaging polarimetry 

The polarimetric data were combined in order to generate a map of the linear polarisation of 
G9. 62+0. 19. We chose the standard approach for obtaining the polarimetric quantities, using the four 
intensities measured with the WoUaston prism (see Sect. l2.2l ). The Stokes parameters can be defined 
as: 

Q = 1(0°) - 1(90°) , (la) 

C/ = 7(45°) -7(135°), (lb) 
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/ = (7(0°) + 7(90°) + 7(45°) + 7(135°))/2 . (1 c) 

The degree of linear polarisation (pun) and the respective polarisation angle (©un) can then be derived 
by using: 

Piin = 0iin = 2 ^ctan — . (2) 

In Fig. IC.8l we show the resulting polarisation map superimposed on the Stokes-7 image (see equa- 
tion II cl) derived from the very same data set. We mention that the degree of linear polarisation as 
derived from Eq. ^ poses an upper limit for the true value of pn^ because of the non-standard error 
distribution of this quantity (e.g., Simmons & Stewartp985). Thus, before plotting we debiased pii„ 
using the approach given in the appendix to the paper of Wardle & Kronberg (119741) . 
NIR imaging polarimetry of circumstellar matter proved to be a valuable tool for revealing the illumi- 
nating source, even in cases where it is hidden from the direct view due to enhanced extinction (e.g., 
Tamura et al. (19911 Yao et al. l2000l Stecklum et al. 2003). In the north-west of Fig.Ol we find two 
adjacent regions of strongly enhanced linear polarisation, covering the positions of component C (cf. 
Fig. IC.2I) and of the H2 emission feature (cf. Fig. IC.7I) . The debiased polarisation degree partly ex- 
ceeds 50 %, indicative of single-scattering events. At first glance, the arrangement of the polarisation 
vectors implies a common illuminator for both regions which seems to be located south-east of them 
towards the more obscured region of the neighbouring IRDC (see Fig. IC.ll) . The ellipse in Fig. IC.8I 
indicates the Icr confidence region where the illuminator is probably located. It is derived by tracing 
the intersections of lines perpendicular to the polarisation vectors whereby we only took the vectors 
with piin > 20 % into account. The most probable location of the illuminator is the centre of mass 
(CoM) of the intersection points, whereas their distribution yields an error estimate. A 2-dimensional 
Gaussian is fitted to the distribution of intersection points, the computed cr parameters, which define 
the fitted Gaussian, are used as axes of the plotted ellipse (Weintraub & Kastner [19931 Stecklum et al. 
12003 1) . The ellipse is quite eccentric because the available polarisation vectors reside in only 2 of the 
4 image quadrants. When looking carefully at the polarisation pattern one sees that the polarisation 
vectors related to the H2 feature are more or less oriented north-south. Hence, when judging by eye, 
one would expect the centre of mass of all intersection points to be located more to the south-east of 
the plotted ellipse. An explanation for the different location in the plot (Fig. lC.8l) is that we can use 
some 180 polarisation vectors near region C for the CoM estimation. But only some 30 vectors from 
the H2 feature region can be taken into account that fulfill the conditions of sufficiently high degree of 
polarisation and acceptable signal-to-noise ratio in the intensity image. Thus, the polarisation vectors 
around component C will contribute far more intersection points than the vectors of the H2 feature and 
hence will simply "preponderate" the CoM estimation. Is this a hint for several illuminating sources? 
We mention that the error Ap is still around 8 % in the high polarised regions. For ratios p/Ap > 3 the 
typical error for the position angle can be approximated by 28.65° x Ap/p (Serkowski 119621) . This 
results in an error of 4° - 9° in our case. But this is of course a statistical error. It could not explain the 
coherent rotation of a whole polarisation pattern in one part of the image (i.e., the H2 feature region). 
The bottom line of these considerations is: The observed polarisation map does not fully correspond 
to one coherent circularly symmetric pattern as expected in the case of only one illuminating object, 
thus indicating the possibility of a second illuminator. 

The radio object E (Fig lC.2l) is situated inside the ellipse in Fig. lC.8l and is an illuminator candidate, 
at least for the scattered light near component C. However, the position of E does not really fit with 
respect to the polarisation pattern of the H2 feature. We mention that the hypercompact radio compo- 
nents G, H, and I (Fig. lC.2l ). whose nature is not yet clarified, are located very near the Icr confidence 



13 



region for the illuminator position and their respective positions are in better agreement with the po- 
larisation pattern of the H2 feature. Thus, one of these sources might be responsible at least for the 
illumination of the H2 feature. Since we do not see clearly corresponding NIR counterparts for these 
objects, they are probably deeply embedded. Should one of these objects turn out to be an illumi- 
nator, this would have implications for the distribution of its circumstellar material. A pronounced 
asymmetry of the ambient matter distribution would be necessary to explain how the object is able to 
illuminate the highly polarised regions in the K band, while we do not detect any K band emission 
from it in our line of sight. 

The K narrow-band and polarisation data of component C prove their importance also regarding the 
general astrometry of the region. In Fig. IC.2I we see that for component C the overlaid VLA con- 
tours of the UCHii region are shifted from the NIR counterpart. One could suspect an error in the 
overall astrometry, if one interprets this K band emission either to arise from Br 7 emission of the 
ionised gas and/or as thermal emission from hot dust, which usually is well mixed with the ionised 
gas of the UCHii region. Our data show that the Br 7 fraction of the K band signal of component C is 
marginal and that the dominant part of this radiation probably consists of scattered light. Under these 
circumstances, there is no reason that the cm continuum contours match the near-infrared emission. 

4. Discussion 

4.1. What is the nature of the IR objects in the HMC region? 

In the previous sections we could establish the newly discovered object F4 as the actual HMC infrared 
counterpart by means of thermal and mid-infrared observations and a careful astrometry. In this sec- 
tion we discuss the consequences that arise from this finding. Some more theoretical consideration 
might help to support our view. 

4.1 .1 . Radiative transfer computations 

To understand the properties and unusual features of the G9. 62-1-0. 19-F hot core we performed some 
basic continuum radiative transfer computations, utilising two selfconsistent radiative transfer codes 
(Manske & Henning |l9981 Wolf et al. 1999|. For this paper we will only consider spherically sym- 
metric models, in accordance with the standard description of hot molecular cores (e.g., Osorio et al. 
11999 ). just to demonstrate the limitations of this concept in the case of our HMC. As input parameters 
we used data from Hofner et al. (.1996J . They had limited the spectral type of the central source to 
be BO. 5 - Bl. The luminosity was estimated to be ~ 1.8 x 10"^ Lq (see also Cesaroni et al. |199 4i 
The mass is not well constrained by other authors, but ranges between 55 and 160 Mq (Hofner et al. 
119961 ). Several power laws for the radial density gradient were applied (p ~ r^" with a = 0.0 . . 2.0). 
As grain material a mixture of silicates (Dorschner et al. 1995 ) and carbonaceous materials (Preibisch 
et al. 119931) was used. The computed fluxes comprise the range of 10"*^ - 10 '^ Jy for A = 2.2 f^m and 
10~'^-5xlO"' Jyfor/1 = ll.Vjum, respectively, and thus are at or below the detection limit of ISAAC, 
while the situation for SpectroCam-10 at 1 1.7 /im is more relaxed. While parameter combinations 
could be found that more or less fit the observational data points for F4 beyond 10 fim, it was not 
possible to also fit the data up to 5 yum. 

In Fig. lC.9l we show the result of such a radiative transfer computation and compare it to the measured 
fluxes of the objects Fl and F4 mentioned in Sects. ITTl 13. 2L and l3.31 The particular parameter set used 
for Fig. lC.9l was chosen to produce preferably high fluxes in the 2 -5 /im range. However, the synthetic 
SED should still be in accordance with the data points at longer wavelengths. For these points, the 
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best fit was reached with a constant density distribution, an outer radius of 2600 AU, a total luminosity 
of 1.88 xlO'* Lq, and a dust mass of 0.94 Mq inside the model space (see Fig. lC.9l ). The optical depth 
for this particular configuration was Tv = 673. It becomes evident that the F4 emission in the 2-5 
/xm range is still more than an order of magnitude stronger than predicted by the model. We interpret 
this as support for the idea that spherical symmetry is not applicable in the case of G9.62+0.19-F. 
However, a comprehensive two-dimensional modelling of the observed fluxes would be beyond the 
scope of this paper. 

To guide the eye, we added in Fig. IC. 91 a modified blackbody curve to the radiative model SED that 
also fits quite well the data points at longer wavelengths. The effective temperature of such a modified 
blackbody would be 90 K, the wavelength dependence of the opacity is decribed by Km ~ /l~' ° in 
our case. However, the significance of such a temperature estimation is questionable for the object 
class of HMCs. As shown, e. g., by Yorke (119801) and Yorke & Shustov (1981) . the concept of an 
effective temperature and an effective radius is not well defined for very dense dusty envelopes around 
protostellar objects, they can strongly deviate from the blackbody or greybody nature. Furthermore, 
we should note that such an effective temperature might be misleading in the following regard. It 
was shown by means of simple radiative transfer considerations, e. g., by Schreyer et al. (il996i) that 
the mass-averaged dust temperature, i. e. the temperature averaged over density and particle size 
distribution, is always lower (AT = 15.. 30 K) than the theoretical colour temperature derived from 
the SED in the case of dense shells around massive YSOs. Thus, most of the mass in such a shell will 
attain a clearly lower temperature than insinuated by the colour temperature. This is important when 
using a temperature in order to derive dust masses from infrared or mm observations. 

The estimated luminosity of 1.88 x lO'* L© from our simple model is in good agreement with 
previous predictions (Cesaroni et al. 19941. When we also take into account the measured excess 
flux of F4 at shorter wavelengths (in comparison to the radiative transfer models and the modified 
blackbody approximation) we gain another 1000 Lq. Hence, the entire luminosity is approximately 
20 000 Lq. This seems to prove that the HMC in G9.62+0.19-F indeed harbours a young high-mass 
star. However, caution is advisable, as we see in the case of the Orion Hot Core, where on a much 
finer scale than in our case several infrared sources can be distinguished (Dougados et al J1993l) . and 
the discussion about the energy budget of the individual sources has not settled yet (e.g., Gezari et al. 



4.1 .2. Hints from the polarimetry 

The polarimetric data can provide additional information on the geometry of the observed regions. In 
a scenario where an outflow has opened the dense molecular shell of a molecular hot core, infrared 
radiation can escape through a channel of reduced density. The light will be scattered at the cavity 
walls ("dust mirroring") which represent an increased column density of scatterers in the optically thin 
cavity environment (Yao et al. 1 19971 120001 ). These single scattering events can lead to high degrees of 
linear polarisation (e.g., Fischer et al. 119961 ). Hence, we would expect to perceive a clear increase of 
NIR polarisation. 

In contrast, the measured polarisation for object Fl is relatively low (10 + 5) %, especially when 
it is compared to polarisation degrees in its immediate vicinity.^ There, the polarisation increases up 
to (33 + 10) % in the small region of diffuse emission to the southwest of F2, but not including F2 



^ It is difficult to confimi this claim just by looking at Fig. IC.8l This kind of plot is well-suited to visualise the 
polarisation pattern for extended emission but is often misleading around compact sources with steep intensity 
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or Fl (best to be seen in the 2.16 fxm image of Fig. IC.4I at the offset position (-0.3, -0.8). Based 
on the low-resolution NIR imaging of Testi et al. (1998), we suggested in our recent outflow paper 
(Hofner et al. 120011) that object Fl itself represents infrared emission arising from the HMC region 
which might be cleared in parts by the molecular outflow. Taking our high-resolution infrared data, 
the improved astrometry and the results of the polarimetry into account, we no longer believe that this 
model can satisfactorily explain the properties of object Fl . 

However, this scenario is certainly a good interpretation for the diff"use emission behind the fore- 
ground star. It appears unlikely that this difi'use emission is just the result of the reflection of light 
coming from the foreground star F2. In such a case of backscattering the reflected light should more 
closely resemble the colour behaviour of the illuminating star. Moreover, there should be a tendency 
for a slighly more bluish SED of the backscattered light due to the decreasing scattering efficiency for 
increasing wavelengths. This is not observed. Instead, the diff"use emission is clearly redder than the 
star and seems to be associated with object F4 whose emission strongly increases in the L' and nbJM 
band, while the star F2 strongly fades beyond 2.2 /um. We therefore think that this star F2 is clearly 
detached from the actual star-forming region. 

We should re-emphasise that both, the weak cm continuum emission and the NH3(5,5) emission as- 
sociated with the HMC component F, peak very close (< 0'.'4, cf. Table EJ to F4 (see also Figs. lC.2l 
and lC.3l) which suggests an intrinsic entanglement between both components. Thus, we interpret F4 
as arising from the inner parts of the HMC outflow cavity. 

One detail is not yet clarified: Is the L' and nbJVI band signal of F4 caused by direct thermal dust 
emission of hot grains very near to the embedded power source of the HMC? Or can a considerable 
fraction of the overall emission be attributed to reprocessed light scattered on the cavity walls? In the 
latter case one has to assume the existence of larger dust grains in the HMC envelope to ensure a 
reasonable high scattering efficiency"' for wavelengths of 3 - 5 fxm. 

These ideas about larger grains are more than a wild speculation. Grain growth is expected to occur 
in YSOs and their near vicinity, in particular, in circumstellar disks (e.g., Beckwith & Sargent 1991] 
Beckwith, Henning, & Nakagawa 2000, D'Alessio, Calvet, & Hartmann l2001l) . Recent observations 
(e.g., McCabe et al. 120031 Shuping et al. 120031) seem to indicate - cum grano salis - the existence of a 
considerable amount of micron-sized particles in disks around low-mass YSOs. Systematic investi- 
gations for the case of massive YSOs seem to be more scarce, especially regarding HMCs. However, 
for the best-known region of ongoing massive star formation, Orion, polarimetric observations at 2.2 
and 3.8 /im point to the existence of larger grains in the Orion Molecular Cloud in general (Rouan 
& Leger fT^MI) and in particular in the BN/KL region containing the Orion Hot Core (Minchin et al. 
[TWn Dougados et al. 1993). 

Thus, a future task will be to conduct high-resolution polarimetric measurements in the K, L', and 
nb_M band. Adaptive optics systems like CONICA/NAOS at the VLT seem to be most suitable to 
ensure high resolution (0 '.'1) combined with high sensitivity - characteristics that are absolutely nec- 
essary in order to deliver meaningful polarimetric results in the case of G9. 62-1-0. 19-F. 



gradients. In these cases the polarisation has to be derived by aperture photometry centred on the compact 
sources in all the polarimetric frames and then applying the Eqs. Q and Q- 

^^ Usually, the scattering efficiency is strongly diminished at wavelengths beyond « 2 - 3 fim when considering 
only dust with the standard MRN size distribution (Mathis et al. 119771 typical for the interstellar medium, i.e., 
the overwhelming majority of the dust particles is sub-micrometer-sized with a typical diameter of « 0.1 jum. 
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4.1 .3. Ordinary stars within tine HIVIC vicinity? 

We have pointed out already, that object F2 very near the HMC is probably a foreground star. Is it 
also possible that object Fl is a normal star whose photospheric emission we see? We have several 
reasons to decline this possibility. 

We start with the following estimate: We use the well-known correlation between the colours of a star 
and its colour excess due to additional reddening: 

H-K = (H-K)o + £h-k. (3) 

Herein, (H - K)q is the unreddened intrinsic colour of the star and the colour excess £'h-k can be 
expressed as 

£h-k=Ah-Ak = (^-1)Ak. (4) 

According to Mathis CI 990"), the AhMk ratio for the standard interstellar reddening law (Ry = 3.1) is 
around 5/3. We can use the intrinsic stellar colours (e.g., Wegner \l 9941 for O and B stars, Bessell & 
Brett 1198^ for stars of type A - M) and our measured value H - K = 4.4 mag to derive the K band 
extinction Ar. O stars of all luminosity classes have {H - K)q values between -0.08 mag and -0.02 
mag. In combining Q and Q we find that Ak is always larger than 6.6 mag for our case. Together 
with the distance module of 13.8 (assuming d = 5.1 kpc) and the measured K band magnitude of 
13.6 mag (see Table El), the absolute K magnitude Mk would be -6.8 mag. Moreover, if we take 
into account that all O stars have (V - K)q values between -0.7 mag and -0.9 mag, the absolute 
visible magnitude My would be brighter than -7.5 mag! O and B stars do not feature such "bright" 
My values, only very extreme supergiants (luminosity class la-0) would exhibit such a characteristic. 
Furthermore, if we considered late-type giants we would also have to deal with extreme luminosity 
classes which, in addition, would cause a severe dilemma to explain the existence of a very old star 
almost in the centre of a young star-forming region. 

Striking evidence against photospheric emission comes from the K band medium-resolution spec- 
trum (Fig. lCTOb which was taken with SOFI at the 3.5-m NTT on La Silla, Chile (Testi l200Tl) . The 
slit was centred on the object Fl, but the spectrum probably covers also minor contributions from 
the foreground star F2. Around 1.87 //m the transmission of the atmosphere strongly decreases due 
to strong H2O absorption, herewith separating the near-infrared H and K band. The strong feature 
around 1.87 yum probably arises from the fact that this atmospheric imprint was not totally canceled 
out in the calibration. In the range between 1.95 and 2.40 //m we see a spectrum without prominent 
features. Especially late-type stars with their relatively cool atmospheres should show many metal- 
lic lines and even molecular lines. A characteristic feature for such stars are more or less strong CO 
absorption features between 2.29 yum and 2.40 yum (Ramirez et al. 119971 Bieging et al. 2002), but we 
only see an almost featureless spectrum monotonically rising towards the K band edge, which is an 
imprint of continuum dust emission. 

These points speak clearly against the "naked star" interpretation for the object Fl. They rather 
indicate the presence of an embedded object within a dusty shell. In a way, this resembles the situation 
of the Orion BN/KL region where several IR objects are located in the vicinity of the Orion Hot Core. 
It might be tempting to see our object Fl as an analogue to the BN object with regard to the general 
appearance and the dusty shell nature of both objects. However, this analogy does not hold. Fig. lC.9l 
indicates that Fl reaches its maximum emission at a wavelength around 5 micron and then drops 
again and hence is not detected at 11.7 yum anymore. This suggests an object of lower luminosity (a 
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few tens L©, if really located at d = 5.7 kpc). BN is 70 - 100 times more luminous than Fl (1500 - 
2500 Lq, Dougados et al. 1993, Gezari et al. 1998 ). it has probably a far more massive reservoir of 
dust around and is of course still strong at 10 and 20 //m. 



4.2. Comparison with other HMCs 

To date, the number of firmly established HMC sources is still quite small. Apart from the well- 
investigated Orion BN-KL region, only a few HMCs have been studied by means of high-resolution 
infrared observations (cf. Stecklum et al. .2001. .2002. De Buizer et al. I2002L 120031 Pascucci et al. 

mm . 

In this respect, it is worth comparing the infrared properties of our region with a similar examination 
of the hot molecular core in W3(H20). Both HMC regions seem to be related to outflow activity. But 
in the latter case no NIR/MIR emission could be detected in the direct vicinity of the HMC (see 
Stecklum et al. 200^. One explanation could be that the outflow is mainly in the plane of the sky, so 
that we cannot benefit from an outflow-related clearing effect. 

Another possibility might be that the potential outflow material had not yet have time enough to 
sufficiently penetrate the very dense shell of the inner HMC region as it is instead apparently the case 
for our HMC component F. In that case, we would rate the 09.62+0. 19 hot molecular core to be in a 
more evolved state than the W3 one. 

Support for this argumentation comes from recent 3-mm and 1.3-mm interferometric line observa- 
tions of 09. 62-1-0. 19 with OVRO (Liu 2003). This author discusses the use of several molecules as 
chemical clocks. First, he uses H2S to re-confirm the age differentiation among the components D, 
E, and F, since H2S is eventually transformed into other species after it has been liberated from grain 
mantles into the warm molecular core environment. Among the three cores, component F shows the 
highest H2S abundance and would thus be younger than E and D which both have already developed 
an ultracompact Hii region. This finding is in accordance with our previous results of the NH3/CH3CN 
abundance comparison (Hofner et al. 119961 ). Furthermore, Liu (120031 ) compares the CH3OH and 
C2H5CN emission at 3 mm. While the observed methanol emission shows extended structure along 
the molecular ridge including D, E and F, ethyl cyanide is primarily detected toward core F. Since re- 
vised chemical HMC models usually produce significant amounts of N-bearing molecules only after 
0-bearing species have reached their peak abundances (e.g., Rodgers & Chamley 2001, 2003), the 
strong C2H5CN emission toward component F suggests that it is in a more advanced stage of HMC 
evolution. 

This fits well to the scenario advocated by Kurtz et al. (120001) that HMCs can evolve first to hy- 
percompact (Oaume et al. 1995, Tieftrunk et al. 119971 Kurtz EUnUI) and finally to ultracompact HII 
regions. The detection of weak emission at 1.3 and 3.6 cm, associated with a high emission measure 
of 3 X 10^ pc cm"^, by Testi et al. (120001 ) implies that the HMC component F is on its way to cultivate 
a hypercompact HII region. 

Furthermore, Cesaroni et al. (119981 ) analysed two other well-known hot molecular cores 
(01 0.47 -1-0.03 and O31.41-l-0.31) and found indications that the HMC stage marks the phase of trans- 
formation from spherically symmetric to flattened and non-spherical structures. Also in this regard, 
09. 62-1-0. 19-F proves to be a typical example of an HMC. 
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Table 4. Physical parameters of the radio components D and E 



Physical parameter 


component D 


component E 


extinction law ^ 


TbDQ ^ 


977 


210 




re[K] 


8000 


8000 




A^(H) [cm-2] 


2.54 X 10^4 


1.84x10^4 




5Brr[Wcm-2] 


1.083x10-1'^ 


2.21 X 10-20 




A2.2^m [mag] 


125 


90 


Ry - 3.1, Cf (2.2yum)-0.108 


A 12pm [mag] 


32 


23 


/?v-3.1,CF(12/im)=0.028 


M.itim [mag] 


232 


168 


Ry = 5.0, CF(2.2A;m)-0.125 


A 12pm [mag] 


60 


43 


7?v - 5.0, CF(12A<m) =0.032 



The brightness temperatures were derived taking the emission in the inner 0'.'55 of the components into 
account which corresponds to our Kg band resolution element (i.e. seeing). 

see Sect. l4.3.^ for details 



4.3. Why can't we see the D and E components in the infrared? 

Radio component F which is only a weak cm continuum source nevertheless shows infrared features 
in its immediate vicinity. Hence, at first glance it might be intriguing that the stronger radio sources 
D and E, which both are more massive and energetic than component F (Hofner et al. 1996), do not 
show an infrared counterpart. Without a further analysis one might come to a similar conclusion as De 
Buizer et al. (120001) who, for lack of an exact astrometry, just shifted their mid-infrared source to the 
position of the radio D component. To enlighten these circumstances, we reconsider radio data from 
previous papers about G9.62+0.19 in order to derive the expected Brackett y fluxes from as well as 
the extinction towards D and E. 

4.3.1 . The expected Bry fluxes of D and E 

The method to compute the expected flux in IR recombination lines on the basis of cm continuum data 
has been explained extensively by Watson et al. (■1997J and was also applied in some of our previous 
investigations of other massive star-forming regions (e.g., Feldt et al. [Tp98 , 1999). We refer to these 
papers for details and start with their formula for the emission measure (EM), which can basically 
be derived from the approximation of Altenholf et al. (I1960I ) for the optical depth of cm free-free 
emission: 

EM = 4.72 «-' T^-^^ v-^ ln( ^—) cm-^ . (5) 

^ e ~ ^ b 

Herein, T^ and T^ are the electron temperature and the brightness temperature, respectively, both 
in Kelvin. The parameter v is the frequency (in Hz) at which T^ was derived. Finally, a{v, T) is a 
correction factor of the order of 1, which is tabulated in IVlezger & Henderson (119671 ) for various 
frequencies and electron temperatures. ^^ 



1 1 Beckert et al. ( 2000 ) demonstrate that this approximation (Altenhoff + Mezger&Henderson correction factor) 
is in very good agreement with non-approximative numerical solutions up to frequencies of ~ 100 GHz. 
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By adopting the approach of Osterbrock (119891) we can use the emission measure to compute the 
expected Br y flux under the assumption that no extinction occurs: 

SBry = 0.9 hVBrr ^r ^ ^^ " (6) 

The factor 0.9 reflects the assumption that a UCHii region consists of ~ 10 % of singly ionised Helium, 
which contributes electrons and ions to the cm free-free emission but only electrons to the Br y emis- 
sion. For a^^ , which comprises the transition coefficients and the level populations via a temperature 
dependency, Hummer & Storey (119871) give: 

a^ff^ = 6.48 X 10"" T;^-^^ . (7) 

For our computations we derived the brightness temperatures of components D and E from the 1.3 
cm map of Testi et al. (2000).^^ The electron temperatures of both regions were estimated by Hofner 
et al. ([1996), using a thermal bremsstrahlung model to fit the cm continuum observations conducted 
up to that point. We list all the parameters in TableHJ 

The expected Bry flux is 1.083 x 10"^'' W cm'^ for component D and 2.21 x lO'^" W cm'^ for 
component E for the theoretical limit of zero extinction. These values correspond to 5.75 mag and 
7.48 mag in terms of K band magnitudes. If we use the ISAAC camera in combination with the Bry 
filter, we would detect such fluxes with a SNR of 20 within less than 1 s of exposure time. 

4.3.2. The extinction towarcJs D ancJ E 

To evaluate these numbers one has to estimate the extinction towards these regions. Hofner et al. 
(119961) derived H2 column densities, based on interferometric C'^O(l-O) observations, for the com- 
ponents D, E, and F They report values of A'^(H2) = 12.7 xlO^^ cm~^ for component D and A'^(H2) = 
9.2 xlO^^ cm"^ for component E. To transform these column densities into extinction values we start 
with the empirical relation that links the total hydrogen column density A^(H) to the optical colour 
excess ^B-i/ (Rvter 1996) : 

A^(H) = 6.83 X 10^^ Eb-v cm'^ mag-^ . (8) 

In order to substitute ^b-Vj we use the ratio Ry of absolute to relative extinction: 

Ry=AylEB^y. (9) 

In combining Eqs. (JHl) and (0) we get an expression for the visible extinction Ay, based on column 
densities. Finally, we can introduce a conversion factor CF{A) that represents the conversion from 
the visible extinction (Ay) to extinction values at other wavelengths (A^). We arrive at the following 
expression: 

^^ = ... tT\ r CFiX) . (10) 

6.83 X 10^' cm 2mag^' 

Mathis (119901) gives two sets of conversion factors, corresponding to two possible extinction laws: 
normal "difi'use interstellar dust" (Ry = 3.1) or "outer-cloud dust" (Ry = 5.0). For both possibilities, 
we compute extinction values at 2.2 yum and 12 //m by using A^(H) = 2 x A^(H2). The results are 



^^ The 3.6 cm map of Testi et al. ( 2000 1 has a higher SNR, but Hofner et al. ( I1996t found that the behaviour of 
the components D and E already deviates from the optically thin regime at this wavelength. 
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listed in Table HJ In any case, the extinction is > 90 mag at 2.2 jum and still > 23 mag at 12 /um. 
Such a large extinction would attenuate the fluxes down to an undetectable value. But one important 
restriction applies to this "naive" interpretation of the extinction data. The C^^O line is optically thin; 
hence, it traces the column density of the whole molecular cloud and not just the column density up 
to the UCHii regions. Theoretically, the UCHii regions could even be located in front of the molecular 
clouds, which is obviously not the case. 

4.3.3. Relative positions of UCHii regions and tine molecular gas 

One approach to clarify the situation is to compare the LSR velocities of radio recombination lines 
(RRLs) produced in the UCHii region with the LSR velocities of molecular lines seen in absorption 
against the radio continuum of the UCHii region. This idea has already been pursued by Downes et 
al. (11980 ) and was recently revived by Araya et al. (2002) who could show that in almost all of their 
UCHii region targets the ionised gas and the main absorbing molecular component have the same 
velocity within the uncertainties of the measurements. This demonstrates that a close association 
between the two classes of objects is very common. 

1.) As far as we know, no radio recombination line studies exist for G9. 62+0. 19 that would provide 
high spatial resolution. Therefore, we go back to the seminal work of Downes et al. (119801) who 
conducted a combined survey of the HllOor RRL and H2CO (lio - In) absorption line with the 
Efl^elsberg 100 m telescope at 6 cm. Their observations for G9. 62+0. 19 were centred on component 
B and had a half-power beamwidth of ~ 2.6 arcmin. The reported LSR velocities are (3.0 + 5.0) km 
s~' for the RRL and (2.0 + 0.5) km s~' for the H2CO absorption line.'^ The velocities are almost 
equal which speaks for a close association also in the case of G9.62+0.19. 

2.) The fact that we actually see absorption implies that a noticeable fraction of the molecular gas is 
located in front of the Hii regions. 

3.) The measured H2CO line temperature for 09. 62+0. 19 is smaller than the continuum temperature 
of the ionised sources (Downes et al. 1980). An opposite finding would imply that most of the H2CO 
absorption occurs just by absorbing photons of the Cosmic Microwave Background and would 
therefore not be related to the Hii regions at all. 

4.) In our HCO+ data of O9.62+0.19 (Hofner et al. EUUll) . which we took with the IRAM 30 m 
telescope and with the Plateau de Bure (PdB) interferometer, we clearly see a strong absorption dip 
in this molecular line at the same LSR velocity as the RRLs and the H2CO absorption occur. An 
interpretation, we mentioned in the paper, is self-absorption due to a high optical depth. However, 
this feature can also be understood as absorption against the 89 OHz continuum of the (ultra)compact 
Hii regions. Probably, both eff'ects contribute to the observed line absorption. Note, that the PdB 
data provide sufficient spatial resolution to distinguish between the components D and E and that 
component B is not dominating in these data. 

5.) Regarding the morphology of the components D and E we do not see features which could explain 
these UCHii regions within the scope of blister/champagne flow models (Tenorio-Tagle I TW^ Yorke 
119861) . which would automatically come along with the location of the UCHii regions very near the 
border of the embedding molecular cloud. 



^^ The RRL velocity is confirmed by Lockman (1989) who gives vlsr - (4.1 ±0.8) km s~^ for the H85q', H87a 
and H88a lines towards G9.62+0.19, based on observations with the NRAO 43 m telescope with a FWHP of 3 
arcmin. 
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These five points support tlie idea tliat both cm continuum components (D and E) are clearly em- 
bedded within their parental molecular cloud, although we cannot exactly estimate their relative posi- 
tions. We cite here recent investigations by Kim & Koo (2003), based on molecular line observations 
towards 16 UCHII regions. They indicate that the formation of massive stars, in general, takes not 
place near the surface of molecular clouds but in their interiors. In this regard, one has to take into 
account that already the positioning of the UCHii regions D and E behind 20 % of the total amount 
of molecular gas would be sufficient to explain our non-detections. On the basis of these considera- 
tions we conclude that the invisibility of the components D and E in the near- and mid-infrared is a 
reasonable finding. 

5. Conclusions 

We have performed a comprehensive set of infrared observations for the G9.62-i-0.19 star-forming 
region. Narrow-band and broad-band observations from 1 to 19 //m provided new information about 
this complex of young stellar objects. In the following we summarise our results. 
1.) The high resolution of our data revealed the detailed structure of the Hot Molecular Core region 
G9. 62-1-0. 19-F. The K band emission found by Testi et al. (1998) could be disentangled into three 
distinct objects. Our astrometry shows that among these three objects the dominating source Fl is 
not coincident with the peak of the molecular line emission of the hot molecular core. Instead, it is 
displaced by ~ 1.7 arcsecs, which translates to roughly 10000 AU on a linear scale. We estimate that 
the second object F2 is probably a foreground star not intrinsically related to the star-forming region. 
No reliable interpretation can be given for the third very faint object F3 in the HMC region, since it 
is only visible in the Ks and L' band. 

2.) Very near the peak of the HMC molecular line emission and partly blended by the foreground star 
we find faint diffuse emission at 2.16 yum. While the foreground star F2 quickly fades at wavelengths 
> 3//m, this emission strongly increases in the thermal infrared and finally reveals the presence of 
another object F4. We propose that object F4 is directly associated with the HMC. This presents the 
first detection of a hot molecular core at a wavelength as short as 3.8 //m. We know that the HMC 
probably drives a molecular outflow roughly aligned with the line of sight (Hofner et al. 120011). Thus, 
the clearing effect of this outflow might allow us to look deeper into the infrared-emitting regions of 
the core environment. The results of our K band polarimetry support the interpretation that the diffuse 
K band emission arises due to scattering of light in the outflow cone. Finally, our mid-infrared data 
highlight object F4 as the most luminous object within the HMC region. With the SED of F4 the 
range for several parameters of the HMC can be narrowed. 

3.) The combined results of the astrometry and radiative transfer computations make it unlikely that 
the dominant object Fl is a direct trace of the HMC. Nevertheless, we have demonstrated that on the 
basis of the NIR colours and the K band spectrum the most likely explanation is a dense dust shell 
around a deeply embedded object of lower luminosity. This might be a hint for multiplicity of YSOs 
within the HMC region. In a way, G9. 62-1-0. 19-F resembles the appearance of the Orion Hot Core 
region, where likewise the near vicinity of the HMC is populated by a variety of infrared sources. 
4.) We evalute the fact that the UCHIIs D and E are not detected in our infrared data although 
these objects show much stronger radio emission than the HMC region F By reexamining available 
high-resolution radio molecular line data we derive theoretical values for the extinction towards D 
and E and show that our non-detections are in agreement with these predictions. 
5.) Using our K-narrow-band results and the imaging polarimetry we reveal within the whole 
G9. 62-1-0. 19 star-forming complex well-defined regions of enhanced Bry and H2 emission as well 
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as a sector where a large contribution comes from scattered light. This demonstrates the complex 
composition of such high-mass star-formation regions that would escape our views when just using 
simple broad-band imaging. 

6.) Generally, we come to the following conclusions. The infrared appearance of complex star 
forming regions is not exclusively governed by the intrinsic spectral energy distribution of the YSOs 
but often dominated by the distribution and asymmetries of the extincting circumstellar material 
on several scales. The interaction of outflows from YSOs with their surroundings can lead to more 
complicated configurations. In particular, infrared emission associated with HMCs is detectable 
under certain favourable circumstances, as we demonstrated for G9.62+0.19-F. 
While spatial low-resolution data, e.g., from the previous generation of IR satellites, could act as 
a principle guide line for investigating massive star formation, high-resolution observations are 
an essential tool to disentangle the objects in high-mass star-forming regions and to make any 
fundamental progress in that field. In this respect, an accurate astrometry is needed since it can have 
a strong impact on the interpretation of the data. 
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Appendix A: NIR Astrometry 

Special care was taken to establish an accurate astrometry for the VLT NIR data. Two approaches 
were used: 

First, we extracted an I-band image from the ALADIN interactive sky atlas (Bonnarel et al. 120001) 
that contained the region of interest. The image we used had been constructed by digital scanning of 
a SERC-l photo plate, using the MAMA/CAI scanning facility (Guibert [19921) ). Such an image has 
- in comparison to DSSl and DSS2 images - a finer pixel scale (O'.'67/px). Within this (ir.5 x 1 r.5) 
subimage we found three useful stars also included in the ACT reference catalogue (Urban et al. 
[1997 ). We calibrated the I-band image with help of these ACT stars and used this image now as an 
astrometric reference to adjust the VLT-Kg image. This could be done straightforward since many 
stars in the central region of the I-band image were also present in the VLT image. 

The second way was more direct without an intermediate step. By applying the ALADIN tool we 
identified more than ten stars in our VLT image also included in the USNO-A2.0 catalogue (Monet 
et al. 1998)^"*. We then used these catalogued stellar positions to directly calibrate the VLT image. 
Nevertheless, one has to be careful when choosing these stars. Some objects which, due to the lower 
resolution, seem to be star-like on the optical photo plate images and are therefore included in star 
catalogues like USNO-A2.0, turn out to be stars with associated nebulosity or clusters of stars when 
imaged with high resolution. Hence, we controled our selection of stars by examining public available 
digitised images of the region. Only positions which appeared to be clearly associated with one single 
object in both, the optical and the NIR images, have been used for the calibration. 

Both methods gave almost identical results. A comparison of several object positions resulted in 
deviations of < 0'.'25 between the two methods. 

After applying these astrometrical calibrations the "red" compact object Fl is still 1.7 arcsec away 
from the position of radio component F - a discrepancy which we find to be significant regarding the 
high resolution of the data. Furthermore, we should mention that the positional difference between 
our astrometry and the one in Testi et al. ( 1998) is obvious, but not large. On the other hand, there 
is a substantial astrometric disagreement between our newly derived positions and the ones given in 
Walsh et al. (fT9^ for their K band image of the G9.62+0. 19 region. 

Persi et al. (2003) use for their investigations of G9. 62+0. 19 the (preliminary) DENIS point 
source catalogue (Epchtein et al. 1999) as astrometric reference. This choice is problematic since 
this catalogue itself uses the old (pre-HIPPARCOS) Guide Star Catalog GSC 1 (Lasker et al.TP^U^ 
for calibration. As a test, we looked for stars in our VLT images which are included in USN02 
as well as in DENIS. Simple marking of the coordinates given in the catalogues reveals that all 
DENIS positions seem to be systematically shifted x 2.5 arcsec to the north-northwest. From that 
spot check we cannot assert a global problem for the astrometry of the DENIS catalogue. But it 
is obvious that the particular stripe of the DENIS release data which contains the G9. 62+0. 19 re- 
gion suffers from an unprecise astrometry and should not be used for calibrating high-resolution data. 



^^ We used the USNO-A2.0 catalogue because it is more or less directly tied to the HIPPARCOS reference 
frame. In the meantime, we also checked the astrometry for the G9.62+0.19 region according to the 2MASS 
All Sky Survey (Cutri et al. 12003 1 . An offset of around 0.3 arcseconds seems to exist between both catalogues 
for G9.62+0.19. We can and will not evaluate here which catalogue is more accurate. We just mention that 
the 2MASS astrometry places our object F4 even more close (O'.'IS) to the peak of the HMC molecular line 
emission (see Fig |C.4l ) 
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Appendix B: Correct position of tlie compact IVIIR source in G9.62+0.19-F 

De Buizer et al. (120001) had first reported on this MIR feature. They used previous 23 GHz continuum 
data (Cesaroni et al. 1 19941 ) to correlate their MIR data. These radio data did not yet reveal the presence 
of the weak 1 .3 cm radio emission arising from component F, but only showed the peak belonging to 
component D (cf. the overlaid contours in Fig. IC.5l) . They then assumed that this radio continuum peak 
for D should be coincident with the MIR peak. However, more recent investigations could already 
show that the compact mid-IR source does not coincide with radio component D (see Stecklum et al. 
l200n Persi et al. l^)03b . Both groups placed the emission within the HMC region G9.62+0.19-F. The 
astrometry was mainly based on the comparison of centroids of the extended emission (component 
B) from the science data and from the (coarser) 12.1 //m MSX images showing the same field. Still, 
the error of such a method is around 2" and does not allow to decide to which of the objects in the 
HMC region the MIR compact emission is related. Finally, De Buizer et al. (2003) presented new 
11.7 fim data that seem to show that the MIR emission is neither associated with component D nor 
with component F but is displaced from both by roughly 2-3 arcseconds. Their new astrometry is 
based on a correlation between MIR features and a 3.5 cm ATCA interferometric radio map taken 
from Phillips etal. (1998). 

In an attempt to clarify the situation, we observed G9. 62+0. 19 at 8.7 yum (Nl filter), using TIMMI2. 
The smaller wavelength (compared to 1 1 .7 /um) and the larger field of view (compared to SpectroCam- 
10) should facilitate the search for compact reference objects in the field. The result is shown in Fig. 
IC.61 where we overlay the TIMMI2 data on our NIR data. In the north, we found a bright star, that 
is also visible in the Kg, L', and M band, suitable for attaching the TIMMI2 astrometry to the VLT 
astrometry. The small rotation misalignment of the TIMMI2 chip with respect to the sky coordinate 
system could be estimated independently, because during the observating run we also observed other 
regions that showed several sources with 2MASS counterparts in the field of view (Grady et al. 120041 
Linz et al. 120041) . The positional uncertainty is around 0.4 arcsecs in RA and DEC (i.e. 2 pixels in 
each direction). We cannot confirm the astrometric result of De Buizer et al. (120031) . We find, that the 
compact MIR emission is indeed located in the direct vicinity of the HMC. Our astrometry speaks 
against a straight identification of the compact MIR emission with object Fl, but rather suggests that 
the peak of the MIR emission is more towards object F4. In comparison to the 1 1 .7 /im data, the MIR 
emission blob is not pointlike. We therefore deconvolved this part of the image with the standard star 
HD 169916 acting as PSF reference. The result (cf. in Fig. lC.4l) is an elongated structure with a similar 
position angle (55° + 6°) as the F1-F4 object pair (51.5° + 2°). This further strengthens our diagnosis 
that also at 8.7 //m we see the objects F4 and Fl that seem to merge to one object due to the lower 
spatial resolution of the the TIMMI2 data. The relative strengths of the contributions of both objects 
is governed by the very different spectral energy distributions of Fl and F4 (see Fig. lC.9l ). 
We should finally mention that one fraction of the astrometrical disagreement between De Buizer et al. 
(12003 ) and our result disappears when agreeing about the exact radio reference position for the HMC! 
We note that they use the J2000 position 18'^06"^ 14^82 (RA), -20°31'38'.'4 (DEC) as HMC reference 
position (cf. their Table 1 and their Fig. 2). But this corresponds neither to the peak of the NH3(5,5) 
emission nor to the position of the faint 3.6 cm continuum of the HMC (see our Table |21 or consult 
the respective references therein). The difference between the 3.6 cm peak position, which we use as 
HMC position throughout this paper and which is the most accurate position for this hot core to date, 
and the De Buizer values is ca. 0.85 arcsec in RA and ca. 1.0 arcsec in DEC, the total displacement 
is hence around 1.3 arcsec. Taking this into account the disagreement between De Buizer's and our 
estimation still does not vanish completely, but has shrunk considerably. 
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Appendix C: MIR image restoration 

If one has, due to the necessary chopping and nodding in the MIR, the positive and negative beams in 
the final image the usual way to treat the data is to perform a kind of shift and add. But in principle, 
it is possible to deconvolve the image using a matrix multiplication procedure. We refer to Bertero et 
al. (200(3 hereafter BBR) for details of this so-called projected Landweber method. The basic idea 
is that within the final multi-beam image several points at different sky positions contribute to the 
signal contained in one pixel of the detector array. Assuming that the chopping and nodding throws 
are colinear along the vertical orientation of the detector chip one can describe the multi-beam image 
g as result of a multiplication of the actually mapped larger rectangular field / with a simple imaging 
matrix A, populated just with natural numbers (see BBR): 

^ = A/. (C.l) 

The projected Landweber method deals with this inversion problem in order to restore the original 
field /. The following iterative ansatz is used: 

^(k.i)=p^[^(k) + ^(A^^_ATA/«)]. (C.2) 

Hereby P+ is the projection onto the positive subset of iterative solutions, A^ means the transposed 
matrix, and r is a relaxation parameter which is related to the largest singular value of the matrix A. 
As starting point /"^"^ = is chosen. 

Applying this method, some effects can influence the quality of the restoration. This contains the ap- 
pearance of ghost images of the restored objects at distances being multiples of the chopping/ nodding 
throws. Such artefacts are inherent in the method itself and can be understood mathematically. The 
BBR article deals with techniques to reduce/remove several kinds of artifacts that may appear. 
Since we used for our observations a chop-nod pattern different from the one used in the BBR ap- 
proach (chopping perpendicular to nodding vs. chopping parallel to nodding), we had to adapt the 
original version of their program. Due to the loss of colinearity the problem is now really a two- 
dimensional one; thus, one imaging matrix alone is not sufficient. The new imaging equation reads 
now with two imaging matrices A and B: 

g = AfB. (C.3) 

The iteration equation evolves to the new key formula: 

y(k+l) ^ p^ [y(k) ^ ^ ^^T^ gT _ ^T^ ^(k) g gTj] (C.4) 

In principle, there is good reason to apply such a kind of restoration method. First, one can easily 
disentangle the regions in the multi-beam image where positive and negative components overlap 
when imaging crowded fields. In this regard, our approach has already proven to be effective 
for the case of TIMMI2 imaging of NGC2264 IRSl (Schreyer et al. EUU^ . Second, one gains a 
larger field of view, because the reconstructed image covers the areas of both, the on-beams and 
the off-beams. This is of course an advantage when one tries to map extended objects with the 
medium-sized (not to say small) MIR detector arrays available nowadays. In addition, within a 
larger field it might also be easier to find other objects that could serve as astrometric reference points. 
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Fig. C.l. Colour-coded image of the entire G9. 62+0. 19 region taken in the three broad-band NIR 
filters J (blue), H (green), and Ks (red). The large-scale contour lines denote the emission levels 
derived from the 8.28 //m image of the related MSX source. The left-most large contour line indicates 
the position of the close-by Infrared Dark Cloud. 
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Fig. C.2. Cutout of the central region of Fig. lC.ll The overlaid contour lines denote the 3.6 cm emis- 
sion of the ionized regions as measured by Testi et al. (2000). The source annotation follows their 
nomenclature. Note, that also component F features weak 3.6 cm emission - the related contour line 
merges with the seeing disk of the yellow star. The black cross marks the peak of the NH3(5,5) HMC 
emission (Hofner et al. .l994j . 
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Fig. C.3. The same field of view as in Fig. IC.2L but now the overlaid contour lines trace the thermal 
NH3(5,5) emission of the HMC as measured by Hofner et al. (1994,). The ellipse in the lower left 
corner indicates the size of the synthesized VLA beam. The black cross at the reference position 
marks the faint peak of the 3.6-cm HMC emission (Testi et al. 2000) . 
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Fig. C.4. Cutout gallery showing the four objects of the HMC region G9.62+0.19-F. Upper row 
from left to right: H, Ks, L' band. Lower row from left to right: nbJVl, Nl, 1 1.7 /im. The size of the 
cutout box is ca. 5'.'4 x 5'.'4. While the foreground star (F2) quickly fades towards longer wavelengths, 
another object very close to it (F4) gets stronger and finally, dominates the region. We presume that it 
is a potential trace of the HMC whose location is here indicated by the upper contours of the NH3(5,5) 
emission . While all the other images are taken from the plain imaging data, the Nl image at 8.7 jum 
shows the TIMMI2 data after deconvolution with the standard star HD 1699 16. Note that the smaller 
spots in that particular image are deconvolution artefacts that change their position and strength when 
varying the deconvolution parameters. 
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Fig. C.5. Inverse gray-scale representation of the 11.7 /urn SpectroCam- 10 image displaying a smaller 
cutout of the G9. 62+0. 19 region. The weaker but nevertheless clearly visible spot to the east repre- 
sents emission from the region of the HMC component F and arises probably from object F4 (see also 
Fig. lC.4l) . Overlaid are the 1.3 cm continuum VLA contours from Cesaroni et al. ( 1994J. The compact 
radio source to the south-east of the compact MIR emission is the UCHii region G9. 62+0. 19 D. 
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Fig. C.6. The JHKs colour composite, overlaid with the TIMMI2 Nl band data as contours. The bright 
source at the offset coordinates (-16, 14) served as main astrometric reference for the TIMMI2 data. 
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Fig. C.7. This image shows G9. 62+0. 19 in the K narrow-band fikers: Bry (blue), H2 (red), (Br y + 
H2)/2 (green). The Bry emission arises mainly from the compact Hii region G9. 62+0. 19 B. Note that 
the H2 emission blob roughly in the image centre has its own substructure. The diffuse emission to 
the north of it (white) turns out to be mainly scattered light (see Fig. IC.8l) . 
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Fig. C.8. Polarisation map overlaid over an image of Stokes-/ (i.e. the intensity) of G9. 62+0. 19, 
both measured at 2.09 yum. The data were binned to a pixel size of 0'.'6 x 0'.'6 before deriving the 
polarisation. The ellipse is an indicator for the location of the illuminator whose scattered light causes 
the high degree of linear polarisation in some parts of the image (see section U31 for details). The thick 
crosses mark the positions of the hypercompact Hii regions E, G, H, and I (cf. Fig. lC.2l ). 
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Fig. C.9. Example for the appearance of the SED of a hot molecular core, calculated with a radiative 
transfer model (solid line). For this graph we used a central heating source with Teff = 22500 K and 
a total luminosity of 1.88 xlO"^ Lq. The total mass contained in the model space was 95 Mq. For 
this particular model we assumed a constant density distribution and an outer radius of 2600 AU. 
No accretion luminosity was added. The drawn-in symbols indicate the fluxes we measure for the 
objects Fl (plus-signs) and F4 (crosses), respectively. The asterisk at 8.7 fim denotes the integrated 
flux from an area covering both Fl and F4. The 2.7 mm flux from Hofner et al. (119961 ) is marked by 
a small square. As comparison, a 90 K modified blackbody is plotted as a dashed line. However, the 
significance of such a temperature is put into perspective in Sect. 14.1.11 
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Fig. C.IO. The K band medium-resolution spectrum of object Fl (in arbitrary units), taken with SOFI 
at the ESO NTT. 



